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We discuss here a model of QCD vacuum in terms of quark antiquark and gluon condensates 
alongwith their fluctuations. The correlation functions of hadronic currents in such a vacuum are 
evaluated to extract hadron properties. The presence of fluctuations of the condensates are em- 
phasized. The structure of vacuum is then generalised to finite temperatures to study correlation 
functions at finite temperatures. Finally, we discuss the vacuum structure at finite densities with 
diquark condensates in a Nambu JonaLasinio type model. The coupled mass gap and supercon- 
ducting gap equations are solved selfconsistently. For certain parameters of the model nontrivial 
solutions of both the gap equations are obtained. The equation of state is also computed. 
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Quantum chromodynamics (QCD) is now accepted to be the theory of strong interaction in terms of quarks and 
gluons, and, at a secondary level, of hadrons. It is a nonabelian gauge theory with SU{2>) as the gauge group. This 
I [ nonabelianess of the interaction leads to two important consequences. Firstly, the interaction become weak at high 
^ ' momentum (above several GeV) transfer processes where perturbation theory is applicable. At low energies and 
00 ' momenta, relevant to most of nuclear physics, the quark gluon coupling strength becomes large and an expansion 
T-H : in powers of this coupling is not useful. The basic difficulty appears to be an understanding of the ground state 
' properties of QCD or its vacuum structure which plays an important role for the related physics . 
, Conceptually, low energy QCD has many common feature with condensed matter physics. The vacuum here appear 
' to be consisting of having fluctuating quarks and gluon fields with average properties being described by condensates 
, of quarks and gluons. The quark condensate (gg), i.e. expectation value of quark scalar densities, plays an important 
' role in the context of strong interactions. It displays the breaking of chiral symmetry. This symmetry is crucial to 
our understanding of nucleons, nuclei and dense matter. The other quantity that charcterise the QCD vacuum is the 
' . gluon condensate. The energy density of vacuum is lowered by presence of electric and magnetic gluon fields in the 
O ' ground state. These condensates were introduced in the context of QCD sumrules and the values are estimated from 
:3 ' the charmonia spectroscopy as (a^TrGG) = {2as/TT^{'E? - B^)) = 0.01206^* and {qq) = -[2?>QMeVf §. 
^ Thus QCD vacuum poses a rich complicated many body problem. As in other many body system ground state 

^ , correlation functions give an insight to the ground state structure. One might be interested in asking questions like 
what happens to a meson pair when placed in such a condensate medium or what happens to interquark interaction as 
a function of their spatial separation. We might remind ourselves that nucleon scattering phase shifts gives information 
regading inter-nucleon forces complementary to that infered from their bound state namely from the properties of 
" " ' deuteron. N-N scattering allows one to study different components of nuclear forces (spin-spin, tensor,- • •) at different 
spatial separation in much more detail than deuteron observables which is an composite effects of all channels. In 
QCD we however do not have free quarks or gluons due to confinement. None the less, one can infer about inter-quark 
interaction as a function of their separation in different channels by studying propagators and correlation functions of 
hadron currents in a more detailed manner than the composite effects reflected in hadron bound state. The correlation 
functions that we shall consider here are spacelike separated correlation functions. To be precise we shall take the 
correlation functions to be defined at equal time so that their separation is purely spatial. Such correlation functions 
have several appealing features. They describe different physics at different spatial separations; they can be calculated 
in some channels phenomenologically from e^e~ — >hadrons data or from r decay experimental data. Moreover they 
can be evaluated in lattice simulations or in some model for vacuum Q . 

Our appraoch here shall be assuming a structure for the vacuum and then examining its consequences regarding 
correlator phenomenology. More precisely, we shall use phenomenological results of correlation functions of hadronic 
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currents [|j to guide us towards a "true" structure of QCD vacuum. We organise this note as follows. In section II, 
we shall discuss a construct for the vacuum in terms of quark and gluon condensates 1^ , and the resulting correlation 
functions |^ . It appears that condensates alone do not give rise to correct phenomenology of correlators unless one 
includes fluctuations of such condensates . Inclusion of the fluctuation fields yields correlation functions consistent 
with the phenomenology of correlation functions. Such a structure of vacuum of zero temeperature is then generalised 
to finite temperature in section III. Here we shall also compute the finite temperature correlators to determine 
temperature dependent hadron properties In section IV, we shall discuss the vacuum structure at high density 
and the QCD vacuum with diquark condensate giving rise to color superconductivity. Finally, we summarize our 
results with some remarks and discussions in section V. 



II. AN ANSATZ FOR THE QCD VACUUM AND CORRELATION FUNCTIONS 

A variational ansatz was considered in Rcf. jsj for the QCD vacuum with an explicit construct involving quark 
antiquark pairs and gluon pairs. The trial ansatz was given as 

\vac)^eMBl-BF){Bl~BG)\0) (1) 
where, the Bogoliubov pair creation operators for the quarks and the gluons are given respectively as 

Bl = J g?(k)t(^ . k)h{k)fj{k)dk, (2a) 

and, 

Bj,^ [ <(k)t.9(k)<(-k)dk. (2b) 



In the above qj,ql^ are two component quark and antiquark annihilation operators respectively. The subscript 
indicates that they annihilate the peruturbative vacuum |0) i.e. q'jlO) = = q^O)- is the gluon annihilation 
operator. The operators satisfy the quantum algebra given in Coulomb gauge as 



and. 



q%(k),q°iiky] ^6^^SrsS{k-k'), (3a) 



i,"(k), a5(k')] = S^" (<5., - k.k./k^) 6ik - k'). (3b) 



Finally, h{k) and ^(k) are two trial functions associated with the quark antiquark condensates and gluon pairs 
respectively. Clearly, a construct as in Eq.(|^)has an obvious parallel to BCS theory of superconductivity. Such a 
structure for vacuum eventually reduces to Bogoliubov transformation for the operators. Then one can calculate the, 
energy functional-the expectation value of the Hamiltonian which is a functional of the condensate functions. Since 
the functions cannot be determined through functional minimisations one can parametrise the condensate functions 
with some trial functions sirnple enough to manipulate numerically as well as reasonable enough to simulate correct 
physical behaviour. In Rcf. ||], the following choices were made (with k = |k|) 

^^"^^^'^)= (exp(i.^P)-I)V. (4a) 



sinh g{k)^A exp{~bk^) (4b) 

In Ref. d the energy density was minimised with respect to the condensate parameters subjected to the constraint 
that the pion decay constant and the gluon condensate value {ag/'K){GG) turns out to be the experimental values 
of 93 MeV and Q.Q\2GeV^ respectively. The result of such a minimisation then leads to instabilty of perturbative 
vacuum to condensate formation beyond a critical coupling of a1 = 0.6. For = 1.28, the charge radius of pion 
comes out correctly. The values of A' and R turns out to be ^' ~ 1 and R ~ 0.96/to. Further some of the baryonic 
properties like charge radius of proton, magnetic moments of proton and neutron turns out to be close to their 
corresponding experimental values. Further, the bag constant- the energy difference between the perturbative and 
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the nonperturbative vacuum turns out to be eo -(140M eV)4which appears to be in general agreement with the 
phenomenological value of this parameter. 

With such a description of the vacuum in terms of condensates let us look at the correlation functions and propa- 
gators in a condensed medium. The equal time propagator is given as |^ 

'9a/3(x)-(i[g„(x),g0(O)]) 

= \ 1^ J t'''' '^^^^^ - 7 • k cos 2h{k)] (5) 

(6) 

Clearly, free massless propagator corresponds to /i(k) ^ limit of above equation and is given as 5*0 = — (i/27r^)(7 • 
x/x^). Different components of the propagator given in Eq.(^) can be analysed Q and it turns out to be qualitatively 
similar to those obtained from other nonperturbative calculations like an instanton liquid model for QCD vacuum 
1^,0. Further, a small x expansion of Eq.(||) yields the propagtor as that one would get in the operator product 
expansion in the vacuum saturation approximation |]lO|. 

We shall next consider the correlation functions of mesonic currents of generic form J(x) = qa{x)Taf}qp(x.), where 
a, (3 are spinor indices and F is a 4x4 matrix in Dirac space (1,75,7^,7^75). The equal time correlation fmiction is 
defined as 

i?(x) = ^{vac\J{x)J{0)\vac) + {vac\J{0)J{x)\vac) (7) 

With the definition of condensate vacuum as in Eq.(^ and the propagator given in Eq.(^, Eq.(j^) reduces to 

R{x) = -Tr[S{x)r'S{-x)r] (8) 

where, F' = 7oF^7o. We shall be normalising the correlation functions with that of free massless correlation function 
which is given as parallel to Eq(||) as 

B.o{x) = -Tr[So{x)r'So{~x)r] (9) 

The ratio R{x) / Rq{x) can then be evaluated for different channels. It turns out that correlation function so obtained 
has similar qualitative behaviour as one can obtain from phenomenology in all chnnels except for the pseudoscalar 
channel 1^. Phenomenologically, in this channel there is a strong attraction with the ratio becoming about 100 at a 
separation of about a fermi The calculated value however turns out to be as low as 1.2 around that value. All 
these results depend very weakly on the functional form of /i(k). 

In view of this outcome, it is obvious that some crucial physics is missing from the model of vacuum considered in 
Eq.(|l]) and has to be supplemented by additional effects. In the present framework this means that quark propagators 
alone do not describe the correlation functions and there has to be contributions from irreducible four point structure 
of the vacuum. This can be thought of as a manifestation of fluctuation of the condensates. Thus we have 

Tqa{-x.)qp{-K)q^{0)qs{0) = Sp-f{:x.)Ssai-^)+ ■ ga(x)g/3(x)g-y(0)g5(0) : 

= 5'/37(x)5fc(-x) + i;/3-y(x)E5a(-z&/a;) (10) 

(11) 

where, we have introduced the composite fields E to include the effects of fluctuation with {vac\T,'S\vac) — but 
(ri|EE|n) ^ 0, where, \ fl)is the "new improved" QCD vacuum including the condensate fluctuations. The correlation 
function then takes the form 

R{x) = - [Tr[S{x)T'S{-x)T] + Tr [S(x))F'E(-a;)F]] (12) 

The structure of E field should be such that it contributes mostly to the pseudoscalar channel and should not affect 
the other channels very much. Such a condition restricts the composite field to be of the form 

Eq/3(x) = /^i(7'7^eyfc(/i''(x) + ^2'5a/30(x) (13) 

where, we have introduced the scalar and vector fileds (p and such that 

{n\^\K)cj^{0)\n)=6'^gv{^) (14a) 
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(l]|0(x)0(O)|n) =5s(x) 
From general considerations, we may write down the functions (x) , 175 (x) as 

1 



27r25 



(14b) 



(15a) 



3s(x) 



[^2-^^1(^42;) - ^J'5Kl{^l6x)] 



(15b) 



Using Eq.s(P,P, p^ , p^ ) one can have the expression for the correlation functions in different channels. Explicitly 
different currents in different channels and the contributions from the propagator and the fluctuation fields are shown 
separately in Table |. We have also included here the currents for the baryons. The resulting correlation functions 
normalised to correlation functions obtained by treating the quarks as massless and noninteracting are plotted in 
Fig.|l|. 
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FIG. 1. The ratio of hadron correlation functions in QCD vacuum to the correlation functions for noninteracting massless 
quarks R{x)/Ro{x) , vs distance x (in fermi). Our results are given by the solid curves. The emperical results determined 
from dispersion analysis of experimantal data in Ref. are shown by long dashed lines. The results from lattice calculations and 
instanton liquid model are denoted by dotted and dot-dashed lines respectively. 
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TABLE I. Meson currents and correlation functions 



CHANNEL 


CURRENT 


CORRELATION FUNCTIONS 

Without fluctuations'' Fluctuat 

(Vector ( 


i?(a;) 
_ Ro{x) 
ion contr 
F^) and 


ibution 
Scalar (F'^)) 


V LUi 




[F{x)f + l^e-^''^' 


F^ = 
F-^ = Stt 


*x«5s(2x) 


Pseudoscalar 






F^' = -487r*a;®3v(2x) 
F'^ = A'K^x'^gs{2x) 


Nucleon 




[F{x)f + ^,^e--'"''F{x) 


F^ = -47r*a;^yv(2a;)F(s) 
= 27r*s^ffs(2a;)F(a;) 


Delta 


eabc [u''{x)C^^,v!'{x)\ u^i^x) 


[F{x)f + ^^e-^-''^'F{x) 


F^ = 
F^ = Stt 


4x«5s(2x)F(x) 



^ F(a;) = [1 + 1^.2/(3.)] where I{x) = /J^ (cosfca; - ^) 

To get the hadron parameters from the correlation functions one relates the correlation function through a dispersion 
relation to a spectral density function. For example in the vector channel the dispersion relation for the correlation 
function reduces to M 



1 

n^^(x) = {n\T{jf,{x)J^{o))\n) = ^ dsR,{s)D{s^'^,x) 



(16) 



with, the function D(ra,x) = {■m/4Tr'^x)Ki{mx) is a propagator of a particle of mass m. The function Ri{x) is the 
normalised spectral density function related to the cross section of e~^e~ annihilation to hadrons. In presence of 
interaction, in general we do not know how to calculate the spectral density function . However we do know their 
qualitative behaviour based on experimental information. To extract hadron parameters one parametrises the spectral 
function and then determine the parameters. A useful parametrisation valid at zero temperature is in terms of a Dirac 
delta function at the pole mass of the hadron accompanied by a step function continuum at higher energy. Specifically, 
for p meson the spectral density function is parametrised by 



Rp{s) = 3A2(5(s - A/2) + 1^0{s - so) 



(17) 



TABLE 11. Fitted Parameters 



CHANNEL 


SOURCE 


M (GeV) 


A 


V^iGeV) 


Vector 


Ours 


0.78 ± 0.005 


(0.42 ± 0.041 GeV)^ 
(0.41 ± 0.02 GeVf 
(0.39 ± 0.02 GeVf 


2.07 ± 0.02 




Lattice 


0.72 ± 0.06 


1.62 ± 0.23 




Instanton 


0.95 ± 0.10 


1.50 ± 0.10 




Phenomenology 


0.78 


(0.409 ± 0.005 GeV)^ 


1.59 ± 0.02 


Pseudoscalar 


Ours 


0.137 ± 0.0001 


(0.475 ± 0.015 GeV)" 
(0.44 ± 0.01 GeYf 
(0.51 ± 0.02 GeYf 
(0.480 GeVf 


2.12 ± 0.083 




Lattice 


0.156 ± 0.01 


< 1.0 




Instanton 


0.142 ± 0.014 


1.36 ± 0.10 




Phenomenology 


0.138 


1.30 ± 0.10 


Nucleon 


Ours 


0.87 ± 0.005 


(0.286 ± 0.041 GeV)'' 


1.91 ± 0.02 




Lattice 


0.95 ± 0.05 


(0.293 ± 0.015 GeV)^ 


< 1.4 




Instanton 


0.96 ± 0.03 


(0.317 ± 0.004 GeV)3 


1.92 ± 0.05 




Sum rule 


1.02 ± 0.12 


(0.337 ± 0.0.014 GeV)^ 


1.5 




Phenomenology 


0.939 


? 


1.44 ± 0.04 


Delta 


Ours 


1.52 ± 0.003 


(0.341 ± 0.041 GeV)^ 


3.10 ± 0.008 




Lattice 


1.43 ± 0.08 


(0.326 ± 0.020 GeV)^ 


3.21 ± 0.34 




Instanton 


1.44 ± 0.07 


(0.321 ± 0.016 GeVf 


1.96 ± 0.10 




Sum rule 


1.37 ± 0.12 


(0.337 ± 0.014 GeVf 


2.1 




Phenomenology 


1.232 


7 


1.96 ± 0.10 
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where, M is the bound state mass, A is the couphng of the current to the bound state and sq is the threshold for 
continuum contribution. The fitted parametrs to the correlation functions obatined by us (solid curve in Fig.|l| are 
given m Table || 

As may be evident from this the results are comparable with that of lattice calculation and instanton liquid 
model calculations Thus to be consistent with the data QCD vacuum must not only have condensates but also 
have their fluctuations. 



III. QCD VACUUM AND CORRELATION FUNCTIONS AT FINITE TEMPERATURE 



As is well known [|ri| the QCD vacuum state changes with temperature. Lattice Monte Carlo simulations suggest 
that chiral symmetry is restored around 150 MeV. Thus it is interesting to look at correlation functions at finite 
temperature in the context of behaviour of hadrons around the chiral phase transition |]l2| , ^ . It may be noted 
that there is little phenomenological information in this regime but there are several theoretical studies 
using operator product expansion (OPE) and sum rule methods as well as using instanton liquid model for QCD 
ground state In particular, we shall generalise the method considerd in the previous section to include 

temperature effects. The zero temperature vacuum defined in Eq.(|l|) can be generalised to finite temperature using 
the methodology of thermofield dynamics. Here the thermal average of an operator is obtained as an expectation value 
of the operator over the thermal vacuum ||20|| . The thermal vacuum is obtained from the zero temperature vacuum by 
a thermal Bogoliubov transformation in an extended Hilbert space involving extra field operators (thermal doubling 
of operators) ||2^ . Explicitly, the thermal vacuum is given as 



= exp (^J cik0(k,/3)(g/t(k)g,t(_k) +q,(k)9/-k)) - h.c}j \vac) 



(18) 



where, the underlined oprators are the operators corresponding to extra of Hilbert space. Further, the q'jS are the 
opeartors refers to the fact that they are the quasi paricle operators corresponding to basis defined by the vacuum of 
Eq.(|l]) i.e. qi\vac) = = q\\vac) and finally 9{k,f3) is the function for the thermal Bogoliubov transformation and is 
related to the number density function given as 



sin2 0(k,/3) 



exp(/3e) + 1 



(19) 



where, e(k) is the single particle energy given as e(k) — -^/k^ + m{k)'^. In the presence of condensate the dynamical 
mass is given as m{k) — fctan2/i(fc) As before the equal time propagator can be calculated including the 
temperature effects as 



(vac,f3\ - [ga(x), (7/3(0)] \vac, (3) 



1 1 



2 (27r)3 



^ cos26l(k, (3) [sin 2/i(k) - 7 • k cos 2/i(k)] 



(20) 
(21) 



As in Section III, we shall take a gaussian ansatz for the condensate function sin2/i(k) = cxp(— i?(r)^k^/2, with 
the condensate scale parameter i?, now being temperature dependent. In order to determine R{T) or equivalently the 
ratio S{T) = R{T = 0)/R{T), we first evaluate our expression of the order parameter (the condensate value) at finite 
temperature. In terms of the dimensionless variable rj = Rk, this is given as 



{qq)T=o 



(22) 



where sit? 6{z,ri) — 1/ (exp(ze(?7)) + 1)), with z — [3/R{T) and e(?7) = ??/ cos2ft(r/). 

We can obtain S(T) = R{T = 0)/R{T) if we know the temperature dependence of the order parameter on the left 
hand side of Eq.(|22|). As there are no phenomenological inputs for this, we shall proceed in the following manner 
to take the temperature dependence of the quark condensate. For low temperatures we shall take the results from 
chiral perturbation theory (CHPT) which is expected to be valid at least for small temperatures [E2|. For higher 
temperatures near the critical temperature, lattice simulations seem to yield the universal behaviour with a large 
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correlation length associated with a second order phase transition for two flavor massless QCD. We shall use such a 
critical behaviour to consider the temperature dependence of the order parameter near the critical temperature. 
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FIG. 2. Figure (a) shows quark condensate at finite temperature normalised to that at zero temperature obtained from CHPT 
and Lattice. Figure (b) shows R(0)/R(T) as determined from Fig (a). 

For intermediate regime we shall take a smooth interpolation between the two. The resulting behaviour of the 
temperature dependance of the quark condensate and the ratio R{0)/R{T) determined using Eq.(|2^) are shown in 



Fig.s 

With the temperature dependence of R{T) known as above, the propagator function of eq.(pT|) now completely 
defined. The finite temperature correlation function is now given as, parallel to Eq.(p^) 

i?(x, T) - - [Tr[5(x, r)r'5(-x, T)r] + Tr([E(x))r'S(-x)r]])T. (23) 

The second term in the above corresponds to contributions from the fluctuations at finite temperature that can be 
written interms of the functions gs(x, T), (x, T) as earlier but now being temperature dependent. We do not know 
how to calculate it except for a general property that the effect of the four point structure should decrease with 
temperature. We take here a simple ansatz for the temperature dependence of gv and gs, 



gs,v{^,T) = 



_ ( 



.gs.i/(x,r = 0) 



(24) 



The parameters /i,; are choosen to have the same values as of zero temperature while fitting the mesonic and baryonic 
correlation function. The normalised correlation functions i?(x, T)/i?o(x, T = 0) are plotted in Fig. 4a and Fig 4b 
respectively. 

As expected (on physical grounds) the amplitude of the correlator decreases with increasing temperature. The peak 
of the vector correlator shifts towards the right after T = 0.9Tc. We might remind ourselves that the position of the 
peak of the correlator is inversely proportional to the mass of the particle in the relevant channel |^ . 

To extract hadronic properties at finite temperature, the correlators are parametrised in terms of a spectral density 
function. This is a generalisation of eq. ( p7| ) to finite temperature and is given as , 



3s 
4^ 



tanh 



AT 



(25) 
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p^(s) = A2 5(s-M2) + |^tanh 



4T 



0(3 -So) 



(26) 



1000 




1.5 



1.0 



0.5 



0.0 




0.0 0.5 1.0 1.5 2.0 

X (in fm) 



0.0 0.5 1.0 1.5 2.0 

X (in fm) 



FIG. 3. The ratio of the meson correlation functions at finite temperature to the correlation functions for noninteracting 
massless quarks at zero temperature R(x, T)/Ro(x, T = 0), vs. distance x (in fm). The solid, dashed, dotted and dot-dashed 
lines correspond to temperatures T =0 MeV, T ^130 MeV, T =140 MeV and T =148 MeV respectively. 

Eq.(p6|) corresponds to spectral density function for pseudoscalar channel. The last term in Eq.(p5|) The last term 
in Eq. (Eq) is the scattering term for soft thermal dissociations (mainly through pions) , which exists only at finite 

temperature and can be taken as QS'p ~ 

The mass, threshold and coupling are then extracted and the results are plotted in Fig. ^ for the vector and in 
the pseudoscalar channel. As can be seen from Fig.s^, with increase in temperature, the correlation functions have a 
lower peak indicating lack of correlations with temperature. In the vector channel the mass of the p meson appears 
to decrease for temperatures beyond 120 MeV. The threshold for the continuum also decreases around the same 
temperature. The behaviour with temperature of these quantities is qualitatively similar to that found by Hatsuda 
et al We have also plotted the temperature dependence of the coupling of the bound state to the current 

which decreases with temperature but rather slowly as compared to mass or the threshold for the continuum. The 
temperature dependence of these parameters can be used to calculate the lepton pair production rate from p in the 
context of ultra relativistic heavy ion collision experiments to estimate vector meson mass shift in the medium. 

In the pseudoscalar channel the mass remains almost constant till the critical temperature whereas the thershold 
and the coupling decrase with the temparature. We have found that in the pseudoscalar channel, the contribution 
to the correlation function mostly comes from the fluctuating fields. Further, the temperature behaviour as taken in 
Eq.(p4[) essentially does not shift the position of the peak whereas the magnitude of the correlator decreases. That 
is refiected in the above behaviour of the parameters in the pseudoscalar channel. We may note here that similar 
behaviour of pion mass becoming almost insensitive to temperature below the critical temperature was also observed 
in Ref. ||2^] where correlation functions were calculated in a QCD motivated effective theory namely the Nambu- Jona 
Lasinio model. 
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FIG. 4. The 
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temperature dependence of mass, threshold (So) and coupling for the vector (left panel) and pseudoscalar (right 
Tc =150 MeV. The vertical lines represent the errors obtained while fitting. 



IV. QCD VACUUM AT FINITE DENSITIES 



Let us now go over to explore the ground state structure in high density QCD. UnUke high temperatures, rather 
little is known about QCD at finite baryon densities from first principles like lattice QCD simulation due to technical 
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problems. However, different low energy models of QCD seem to indicate a rich phase structure in this domain. In 
particular the color superconductivity phase at high density has attracted much attention. Although it was known 
for quite sometime [|^ , recent studies [p5|-p8| indicated that the superconducting gap could be as large as 100 MeV 
and that has lead to an extensive literature on this subject in recent past regarding its consequences in neutron stars 
as well as for heavy ion collisions. 

As per BCS theorem, an arbitrarily weak attractive interaction makes the fcrmi sea of quarks unsatble at high 
densities. In deed, the color current current interaction is attractive in the scalar (3) and axial vector channel. 
This can lead to formation of Cooper pairs and associated superconductivity in the color space. To discuss it in a 
nonperturbative manner, we take the trial ansatz over the chiral condensed vacuum of Eq.([^ as Q 



\oj) = exp 



dk. 



(27) 



In the above, i,j are flavor indices, a, 6 are the color indices and r(= ±1/2) is the spin index. We shall consider 
here two flavor and SU(3) color. Here we have also introduced two trial functions /(k) and /i(k) respectively for the 
diquark and diantiquark channel. As may be noted the state constructed in Eq. (|2^) is spin singlet and is antisymmetric 
in color and flavor. 

Next, to include the effect of temperature and density, we obtain the state at finite temperature and densit y \lo{I3 , /i)) 
by a thermal Bogoliubov transformation over the state \uj) using thcrmofield dynamics as in Eq.(O) |p|,po|,pip^] , 



I o;(/3, /i)) = exp l^j qi (k) 1 0_ (k, (3, fi)qj (k) ^ + qj {k)0+ (k, /3, M)g^ {k)dk - h.c. 



(28) 



In the above, the ansatz functions 6±{k, (3, fi), as before, will be related to quark and antiquark distributions. We 



might note here that the trial ansatz given in Eq.(28) actually invloves five functions - h{k), for the quark anti 
quark condensates, /(k) and /i(k) describing respectively the diquark and diantiquark condensates and 0±{k, /3, /i) to 
include the temperature and density effects. All these functions are to be obtained by minimising the thermodynamic 
potential. This will involve an assumption about the effective hamiltonian . For the purpose of illustration we shall 
consider a hamiltonian of Nambu-Jonalasinio type given as 

H = ^^(-^a•V)^ + ^J;:J''^ (29) 

with = q7^r"q. One can then calculate the energy functional e = (a;(/3, /i)|7i|a;(/3, /i)) and the thermodynamic 
potential JF = e — {1/(3)S — fiN, using the Bogoliubov technique. The details are given in Ref. |^ . The free energy is 
a functional of all the five functions h, /, fi,0±. However with the point interaction of the the Hamiltonian of eq.([29|) 
one can determine them. This results in two coupled gap equations to be solved in a self consistent manner and are 
given as 

^9" 1 f dk^^^ ( ^tanh{^) + ^tanhi^)] ^ 1 (30a) 



3 (27r)3 J Vk2 + M2 V^^- 2 ' oj+ 




= 1 



(30b) 



where, uj± = y + and ^± — [E Here, E is the energy of the quasi paricles given as E' = (A;^ + M^), 

and,z/ is the chemical potential in presence of interaction given as = /i — {Ag^ /2>){N / 12), N being the quark number 



density. Eq.(30a) is the mass gap equation in presence of diquark condensates and Eq.(30b) is the superconducting 
gap equa tion which is a relativistic generalisation of BCS gap equation | |30| , |3l[ | . In the limit of no diquark condensates 
Eq.( |30a|) reduces to that obtained in Ref. |^ except for the numerical factors before the integrand. This is due to the 
fact that the approximation the the later case has been a mean field approximation [p2[ unlike the case here where 



approximation lies only with the ansatz for the ground state. Eq.(30a) without the diquark condensate contributions 
is also have the same structure as in Ref. in the limit of the formfactors intoduced in the later case reduces to 
a sharp cutoff in the momentum. Similarly in the limit of chiral condensate going to zero, the super conducting gap 
equation (pOb|) is similar to that obtained in Ref. |M| or Ref. p5|. 
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FIG. 5. Mass gap M as a function of fermi mometum. The dashed line correspond to no diqurak condensates. The solid line 
corresponds to both diquark and quark antiquark condensates structure for the ground state. 

The solution of these equation at zero temperature is shown in Fig.(^ for the mass gap and in Fig.(^) for the 
superconducting gap. For the sake of comparision we have also plotted the mass gap without the diquark condensates. 
The coupling here taken as ~ 56GeV~^ and the cutoffA ~ 0.67GeV. these values are taken so as to give the same 
transition temperature as in ref p5| , pl| . With these couplings, the mass gap at zero temperature and density is about 
490 MeV and the mass gap vanishes at fermi momentum kf = ^/v^ — = 400Mel^. The corresponding critical 
quark number density is about \.l j ]rn? . 
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FIG. 6. Superconducting gap A as a function of fermi mometum. The solid line corresponds to both diquark and quark 
antiquark condensates. The dashed line corresponds to only diquark condensates. 



Presence of diquark condensate does not change these values very much. The diquark condensate increases with 
number density and becomes maximum of about 90 MeV beyond which the effect of the cutoff is felt and it vanishes 
for fermi momentum around 600 MeV. We also plot the equation of state (EOS) in Fig.0. 
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FIG. 7. Pressure as a function of fermi momentum 



While plotting the EOS, we have added the bag constant eq, which is the difference in energy density of the 
perturbative vacuum and the nonperturbative vacuum with condensate at zero temperature and density, to the 
pressure. With the present parameters the bag constant turns out to be — (200MeT^)^. The pressure has a cusp like 
structure and becomes negative at finite densities. The portion of the curve that goes down with kf corresponds to 
the nontrivial solution to the mass gap and the portion that increases with density correspond to zero mass solution. 
The negative pressure indicates mechanical instability and can have the interpretation that uniform nonzero density 
matter will break up in to droplets of finite density in which chiral symmetry is restored surrounded by empty space 
with zero pressure and density [ p5[ . It is tempting to identify the droplets of quark matter with nucleons within which 
density is nozero and {qq) = - a fact reminscent of bag models [p5| . Nothing in the model however says that the 
droplets have quark number three. The equation of state does not change much in presence of diquark condensates. 
This should be expected because the effect of diquark condensate is small in the region where chiral condensate is non 
vanishing. Thus gross structural properties of the neutron stars are not likely to be affected by diquark condensates. 
However, cooling of neutron stars shall be expected to be very much affected by such a gap of about 100 MeV. 



V. SUMMARY AND DISCUSSIONS 



We have looked into the structure of QCD vacuum in a nonperturbative manner with a variational ansatz. This 
has been done for zero temperature as well as finite temperature and densities. The input has been equal time algebra 
for the field operators and the ansatz for the ground state. At zero temperature, the Bogoliubov type pairing ansatz 
involving both quark and gluon condensates becomes energetically favourable beyond a critical coupling. It also gave 
some of the low energy hadronic properties like pion and proton charge raddi, proton, neutron magnetic moments 
and the bag constant to be around their phenomenological values for the coupHng value of 1.28. Thus as = 1-28 
effectively corresponds to the QCD cou plin g constant for the vacuum configuration. With optimized renormalisation 
group equations, it has been seen that [|33| as{Q) does not go to infinity as Q decreases below 300 MeV, but freezes 
to a constant value around unity. Our analysis seems to remind us of a similar situation. 

We next evaluated the correlation functions of hadronic currents in such a condensed medium. However it appears 
that to have quantitative agreement with correlator phenomenology, particularly in the pscudoscalar channel, these 
cannot be described by the propagators alone but must of necessity have the fluctuations of the condensates. This 
may be looked upon as combination of two effects -(i) an effective way of incorporating gluon condensate effects and 
(ii)the existance of explicit four point structure in QCD vacuum. In some ways these fluctuations may be related to 
the "hidden contributions" discussed by Suryak pO| . 

It is worthwhile pointing out that OPE and our approach are based on intrinsically different assumptions. The 
former is an expansion which separates short distance (Wilson coefficients) and long distance (condensates) physics. 
In our method we assume an explicit vacuum structure in terms of quark condensate (two point function) plus an 
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irreducible four point function. Having made an ansatz for the vacuum, we do not make any further approximation 
in the evaluation of the correlators. The approach is phenomenological in the sense that the values of the parameters 
in the four point function 1^ are chosen to reproduce the behaviour of the correlators. As emphasized by Shuryak 0| 
and Schiifer and Shuryak |1§[| OPE is able to quantitatively describe the zero temperature pion correlation functions 
for small x (upto 0.25 fm) but underpredicts it for large x. In our work the agreement is quantitative with 
experimentally deduced mesonic correlation function and lattice results for the whole range. This covers the small 
x values, resonance region and large x domain where the correlator vanishes. We have however more parameters. To 
the extent that the large x behaviour of the pionic correlator depends on gluon condensates in OPE our parametrization 
would imply an effective way of including gluon condensate effects. 

For studying the correlators at finite temperature, we assume that the two point as well as the four point function 
vanish at T = Tc- The parameters in the four point function are assumed constant at their T = values - no 
additional T dependence is given to them. We see that our results are broadly in agreement with those of Hatsuda 
etal [0. All the same, it is not possible to carry out a term by term comparison of our results with OPE results. 
This is because the assumptions are different in the two approaches. We recall that Hatsuda etal attribute 
the decrease in the p mass to contributions coming from four point function. They also find variation in the gluon 
condensates to be less than 5% over the temperature region that is considered. In our work we do not include the 
gluon condensates and the four point structure vanishes as T ^ Tc- This is not the case for gluon condensate {{GG)) 
in lattice calculations or in the dilute pion gas model of Ref. [ p^t . Consequently, we may be tempted to infer that the 
decrease in rho mass in our model is due to "genuine" four point function effects (as in Hatsuda etal |l7| ) and not 
from the "effective" gluon condensate contribution. One, however, has to be very cautious. This is because in the 
present work the parameters in the four point function were kept constant at their T = values. It should be recalled 
that these values were determined so as to correctly describe the behaviour of the correlators (in particular pion) 
at T=0. Hence the parameters do reflect some effective gluon condensate effects. The results at finite temperature 
would certainly be modified if these parameters are given significant T dependence. Our parametrization is such that 
if the parameters decrease with T then the contribution to the correlator will decrease. The crucial question however 
still remains as to the behaviour of (GG) for T < Tc- If it varies very little in this range then our assumption of the 
parameters remaining constant would be reasonable. 

We would like to add here that the present analysis will be valid for temperatures below the critical temperature. 
Above the critical temperature there have been calculations essentially using finite temperature perturbative QCD 
in random phase approximations (RPA) jsj]. However, in the region above Tc, nonperurbative features have been 
known to exist from studies in lattice QCD simulations MM. In view of this, one may have to carry out a hard 



thermal loop calculation where a partial resummation is done 35 1 . Alternatively, one may use other nonperturbative 
approaches such as QCD sum rules at finite temperature or RPA approach in an instanton liquid model for the 
QCD vacuum 

The vacuum structure at finite densities was looked into taking into account the possibility of diquark condensates 
in a Nambu-Jonalasinio model. As earlier, the approximation lies here only in the ansatz. The resulting gap equation 
for the superconducting gap is a relativistic generalisation of the BCS gap equation. The mass gap equation in the 
presence of superconducting gap is a new feature of the present calculation. Because of the point interaction as 
in Hint of Eq.(|2^)we could solve for the gap functions explicitly. In the presence of realistic potentials as in Ref. 
|37|or Ref. one will have to solve an integral equation for the gap functions. Such a calculation is in progress 
|29|. It will be interesting to see how the results for the superconducting gap would be affected in presence of quark 
antiquark condensates as compared to the resummed perturbative QCD calculations at finite densities The 
equation of state here did not change very much. Thus the global properties of neutron stars shall not be affected 
in an appreciable manner. However, a gap of about 100 MeV can have its implications on neutron star cooloing 
pof , magnetic fields of pulsars and their thermal evolution [Q. Future theoretical studies of QCD at finite baryon 
densities may reveal to which extent these newly established features of quark matter considered in effective models 
shall have their correspondence in a more complete treatment of QCD at finite baryonic density and shed light on 
whether one may expect any distinguishing feature in the global properties of neutron/quark stars. 
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